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ABSTRACT 



We present constraints on the evolution of large-scale structure from a catalog of 710,000 
galaxies with Iab < 24 derived from a KPNO 4m CCD imaging survey of a contiguous 4° x 4° 
region. The advantage of using large contiguous surveys for measuring clustering properties on 
even modest angular scales is substantial: the effects of cosmic scatter are strongly suppressed. 
We provide highly accurate measurements of the two-point angular correlation function, u>(6), 
as a function of magnitude on scales up to 1.5°. The amplitude of uj{6) declines by a factor 
of ~ 10 over the range 16 < I < 20 but only by a factor of 2 — 3 over the range 20 < 
/ < 23. For a redshift dependence of the spatial correlation function, £(r), parameterized as 
£(r,z) = (^)~ 7 (1 + z)~( 3+e \ we find r Q = 5.2 ± OAh' 1 Mpc, and e > for I < 20. This is 
in good agreement with the results from local redshift surveys. At / > 20, our best fit values 
shift towards lower r Q and more negative e. A strong covariance between r D and e prevent 
us from rejecting e > even at faint magnitudes but if e > 1, we strongly reject r £ 4/i _1 
Mpc (co- moving). The above expression for £(r, z) and our data give a correlation length of 
r Q (z = 0.5) ~ 3.0±0.4/i -1 Mpc, about a factor of 2 larger than the correlation length at z = 0.5 
derived from the Canada-France Redshift Survey (CFRS; [Le Fevre et al. (1996)| ). The small 
volume sampled by the CFRS and other deep redshift probes, however, make these spatial 
surveys strongly susceptible to cosmic scatter and will tend to bias their derived correlation 
lengths low. Our results are consistent with redshift distributions in which ~ 30 — 50% of 
the galaxies at I = 23 lie at z > 1. The best fit power law slope of the correlation function 
remains independent of / magnitude for / < 22. At fainter limits, there is a suggestive trend 



towards flatter slopes that occurs at fluxes consistent with similar trends seen by Neuschaffer 



& Windhorst (1995) and Campos et al. (1995), The galaxy counts span 11 magnitudes and 
provide an accurate calibration of the galaxy surface density. We find evidence for mild galaxy 
evolution - about 1 mag of brightening or a doubling of the density by I = 23 relative to an 
Q = 1 no evolution model; about 0.5 mag of brightening or a factor of 1.5 increase in surface 
density relative to an open model. Our galaxy counts agree well with those from the HDF 
survey and, thus, argue against a significant inclusion of sub-galactic components in the latter 
census for / < 24. 



Subject headings: large-scale structure, clustering, galaxy evolution, galaxy catalogs 



1. Introduction 

The evolution of large-scale structure (LSS) in the universe probes the abundance and form of dark 
matter, the mean baryon density, the turnover scale in the perturbation power spectrum, and the formation 
processes of galaxies and clusters. Observational constraints on the evolution of LSS provide strong limits on 
structure formation theories because the sensitivity of the mass function to f2 is high and because coherent 
structures on scales > lO/i -1 Mpc are still in the linear regime, precisely where important distinctions can 
be made between competing theories. A study of this evolution is, thus, fundamental. Substantial progress 
has been made in studying LSS evolution from deep redshift surveys. The general consensus from the most 
distant probes flCohen et al. (1996) , Stcidel et al. (1998) , Giavalisco et al. (1997) ) suggests that there has 



been significant evolution in the clustering properties of galaxies since z ~ 3. Precisely just how much 
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evolution has occurred since z ~ 1 is less certain but the data are being accumulated rapidly (Lilly et al. 
(1995) |, |Le Fevre et al. (1995)1 , jConnolly et al. (1996)j |Connolly, Szalay, fc Brunner (1998)| ). 



Two-dimensional surveys offer some key advantages over narrow, deep redshift surveys like those 
above. Specifically, a large area angular survey subtends many galaxy correlation lengths in the transverse 
direction enabling many independent cells to contribute to the signal at a given depth. In addition, the 
power spectra derived from angular surveys are unaffected by redshift distortions, which tend to wash out 
small-scale power. The disadvantage of a 2D survey, of course, is that the depth sampled in a given flux 
range is usually broad and one must therefore deal with projection effects. Observationally, mapping large 
areas of sky is now quite an efficient process due to the availability of large-format CCD arrays and mosaic 
cameras. Ideally, one would like to conduct deep redshift surveys that also cover large, quasi-contiguous 
areas. As an initial step towards this goal, we have surveyed a contiguous 16 square degree area in the 
/-band using the prime focus CCD camera on the Mayall 4m at Kitt Peak National Observatory. The area 
is about an order of magnitude larger than any comparably deep, contiguous survey yet published. Our 
goal is to set accurate limits at z = 1 on the amplitude of structures now seen at z = 0. 

The study of large-scale structure has advanced primarily through the study of the moments of the 
galaxy distribution. These fundamental functions, which must be explained by any viable structure forma- 
tion model, are determined, in part, by the initial conditions at recombination and by subsequent non-linear 
growth. The presence of non-zero high order moments (n > 3) provide constraints on the non-Gaussian 
properties of the galaxy distribution and, hence, provide very strong discrimination between competing 
models. The first two moments, however, are essential for assessing departures from a non-evolving, Pois- 
son distribution. The combination of area and depth provided by this survey (~ 710, 000 galaxies with 
Iab < 24) allows us to determine the first two moments with significantly lower uncertainties than previous 
works. For example, the angular two-point correlation function at faint magnitudes has been measured 
on scales of 10 arcminutes (e.g., Efstathiou et al. (1991), Neuschaffer fc Windhorst (1995)| , |Campos~ 



al. (199(5]], |Lidman fc Peterson (1996)j |Woods fc Fahlman (1997)1 , [Brainerd fc Smail (1997] ) but there is 



substantial scatter in the results at I ^ 21, which has lead to conflicting conclusions about the evolution of 
structure. A significant reason for the scatter is the relatively small ( ;$ 0.25 deg 2 ) contiguous areas used. 
At these scales, one is subject to both cosmic variations in clustering and systematic corrections that are 
comparable in amplitude with the signal being measured. Such effects are negligible for the present survey 
and thus enable a major enhancement to the precision and accuracy of LSS and evolution measurements. 
Furthermore, the survey supports a reliable determination of at least the 3rd and 4th moments of the 
galaxy distribution. 

In this paper we present the /-band galaxy counts and constraints on the evolution of the two- 
point correlation function. An analysis of the higher moments (n = 3, 4) of the galaxy distribution will 



be presented in a separate paper ( Szapudi et al. (1998) ). Section 2 contains a brief description of the 



observations and survey strategy. Section 3 presents results for the /-band galaxy counts. Our results for 
the two-point angular correlation function are presented in Section 4. The interpretation of our results is 
discussed in Section 5. We adopt h = H o /(W0 km sec -1 Mpc -1 ). 



2. Survey Strategy and Observations 



The primary goal of the survey is to use distant clusters of galaxies to explore LSS out to z = 1. To do 
this reliably, the survey had to subtend at least ~ 75/t" 1 Mpc in at least one of the angular directions given 
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that the local correlation length of clusters is 15 — 20/i 1 Mpc ( postman et al. (1992) , Dalton et al. (1994) 



Croft et al. (1997)| ). At z = 1, 75/i -1 Mpc subtends 4.1° for q Q = 0.1 (4.9° for q a = 0.5). Furthermore, 



to measure the high redshift cluster correlation function with accuracies approaching those of low redshift 
studies, the survey needed to contain at least 150 - 200 clusters with z > 0.3. Estimates of the surface 
density of clusters (Gunn et al. (1986] , Postman et al. (1996)| ) are in the range 10 - 15 clusters/deg 2 out 



to z = 1, thus an area of at least 15 square degrees was needed. Additional survey size constraints were 
imposed by the desire to accurately measure the galaxy angular two-point correlation function, iv(9), on 
scales up to ~ 1° at faint magnitudes. This required that the survey subtend at least 3° to minimize the 
effects of the integral constraint on uj(9) on scales 9 < 1°. These requirements, plus the desire to avoid 
severe aliasing effects and unfavorable window functions, led us to choose a 4° x 4° survey geometry. 

The KPNO 4m prime focus CCD camera has a 16' field of view (0.47 arcseconds per pixel); 256 
exposures were, therefore, required to survey the entire field. Each pointing overlapped its adjacent pointing 
by 1 arcminute. Local stability of the zeropoint is provided by comparing the overlap regions of any given 
CCD field with its neighbors. By comparing the photometry of objects in common in any overlap region, 
we find that we can measure the relative photometric offsets between any two images to 0.03 mag. We 
then solve for relative offsets of all images in the survey by a system of linear equations that minimizes 
the local differences between the images and their neighbors. Simulations of this operation show that 
the local error in the zeropoint of any image is ~ 0.016 mag which translates to a limit on the random 
error in uj(9) of ^ 0.001 due to frame-to-frame zeropoint variations. As roughly half of the images were 
obtained in non-photometric conditions, this allows us to still calibrate the survey from the photometric 
set of images. Unfortunately, building the calibration from a set of local differences is poor at preserving 
a constant zeropoint over large scales. To provide for stability of the zeropoint over the extent of the 
survey, we used large-area calibration frames obtained under photometric conditions at the KPNO 0.9m 
telescope by Ralph Shupping. The calibration frames were interspersed throughout the survey area, and 
provide a check against large-scale gradients in the calibration. Based on the quality of this material, 
we estimate that the zeropoint is constant over both 4° dimensions of the survey to ^ 0.04 mag which 
translates to a systematic error in lo(9) £ 0.003 on a 4° scale and proportionally less on smaller scales. 
Indeed, computation of u>(9) in the independent quadrants of the survey yields results which are consistent 
with one another at the above level. 

We selected the field centered at lOh 13m 27.95s +52d 36m 43.5s (J2000) by virtue of its high galactic 
latitude (+51°), low HI column density (2.2x 10 20 cm -2 ), high declination (extended visibility from KPNO), 
low IRAS 100/x cirrus emission, and the absence of many bright stars or nearby rich clusters. In a survey 
of this size, however, it is impossible to avoid all bright (/ < 16) stars. The effective area of the survey is 
14.7 deg 2 after we exclude those pixels in the vicinity of these bright objects. 

Each exposure was 900 seconds in duration yielding images with a 5<r limit at Iab = 24 and sufficient 
depth to detect cluster galaxies 2 magnitudes fainter than the typical unevolved first-ranked elliptical at 
z = 1. This depth was essential — a shallower survey would only be sufficient for detecting the very richest 
z ~ 1 clusters (at Iab = 24, we are able to detect z ~ 1 Abell richness class 1 systems) and would limit 
visibility of structure evolution. 
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3. Galaxy Counts 

Automated object detection and classification were performed using a modified version of the FOCAS 
package (Valdes (1982)). Modifications include the use of a position-dependent point spread function 



(PSF), an essential feature for accurate star/galaxy classification in data obtained with the KPNO PFCCD 
system. Extensive catalog quality assurance was performed to verify object classification accuracy and to 
remove spurious detections near bright objects. A detailed description of the galaxy catalog construction 



process is discussed in Postman et aL (1998) . However, it is important to provide a brief description of 
the star-galaxy classification procedure here. The object classification is done on a frame-by-frame basis. 
First, the PSF in each image is determined using a set of compact, symmetric objects (this is accomplished 
with the FOCAS autopsf routine) . To assure the code has not selected inappropriate objects, each derived 
PSF is visually inspected and a quick comparison is made with the actual image frame to assess whether 
the PSF template is sensible. In about 10 - 15% of the cases, stellar templates had to be selected manually. 
Second, a polynomial representation of the PSF variation as a function of distance, r, from the camera's 
optical center is made. (Specifically, we find the function B(r) = a® + a2r 2 + a^r^ is an accurate model 
of the variation in the KPNO camera prior to the installation of the new prime focus corrector/ ADC in 
June 1997). The optical center is determined directly from the data by mapping the stellar surface density 
(using a position-independent PSF classifier) as a function of CCD position. Failure to use a position- 
dependent PSF results in a paucity of stars near the CCD edges which traces the optical distortions. 
Finally, each detected object is classified as a star or galaxy based on fitting both an unbroadened and 
broadened PSF to the two-dimensional object image. Any object which requires a component broader 
than 1.3 times the nominal PSF at the object's position on the frame is considered to be a galaxy. The 
stellar surface density derived from position-dependent PSF classification shows no variation with distance 
from the optical center. This classification method works well down to about 1=21.5 in these data (based 
on running the same code on simulations). Fainter than this, a statistical approach is needed to assess 
the probability that an object is a galaxy. However, for I > 21.5 galaxies outnumber stars by at least a 
factor of 5 (the galaxy/star ratio is ~ 10 by I = 23.5, based on an extrapolation of the star and galaxy 
counts measured to I = 21.5). Hence, star contamination would not exceed ~ 10 — 15% even if one simply 
classified every object fainter than I = 21.5 as a galaxy. We, however, do use a probability derived from 
the extrapolated galaxy/star ratio to classify objects fainter than / = 21.5. 

Figure |l] and Table 1 show the differential number of galaxies per unit area per magnitude as a 
function of I magnitude. Also shown in Figure |l| are the predictions from the dwarf-dominated, no- 



evolution models of Ferguson fc Babul (1998) and the counts from the Hubble Deep Field (Williams et al. 
(1996)|J]all adjusted to the Cousins I system. The curvature seen in the I— band log N-mag relationship 
is real - the effects of object misclassification and incompleteness are negligible for I < 21 and have been 
accurately modeled and removed at fainter magnitudes as discussed above. The slope of the relationship 
monotonically decreases in a nearly linear fashion from the Euclidean value of 0.60 at I = 15 to 0.27 by 
/ = 23.5. Figure demonstrates how our counts compare with previous determinations. For clarity, we 



have normalized all results by dex(0.6(I — 18)). Our counts are in good agreement with those of Lilly 



et al. ([991], [Smail et al. (1995)| , |Le Fevre et al. (1995)| (CFRS), |Williams et al. (1996)| , and pardner 



et al. ( f996)j Our survey provides a homogeneous calibration of the /—band galaxy counts over an 11 
magnitude range (13 < I < 24), and actually covers more area at the bright end (I < 19) than any of these 
earlier surveys. This is an important step towards removing ambiguities in the galaxy number counts at 
brighter magnitudes. We find that the CFRS counts are systematically high for J < 19 and the Hubble 
Deep Field (HDF) counts are systematically high for / < 21. However, the errors for those surveys are 
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large in those intervals owing to their small angular coverage and the difference is entirely consistent with 
the observed cosmic fluctuations in counts on scales of a few arcminutes and probable small differences 
in isophotal detection limits. Our agreement with the HDF counts in the range 21.7 I 24 is quite 
important because it suggests that both the HDF galaxy catalog and ours appear to be counting the same 
type of objects at least in this magnitude range. Given HST's superior resolution, this would argue that 
down to I = 24 inclusion of a significant population of sub-galactic clumps in the HDF counts is not likely 



flColley et al. (1996) 



We have normalized the no-evolution models to match our counts precisely at I < 16 where the models 
and the data exhibit the same near-Euclidean slope of 0.55. This is the most appropriate normalization 
location since we want to test for departures from the no-evolution assumption at cosmological distances. 
Relative to the no-evolution = 1 model at / = 23 (0.7 £ -^median ~ 1)) our counts are consistent with 
either ~ 1 mag of luminosity evolution or a factor of ~ 2 enhancement in the galaxy density. Relative to 
the £1 = 0.02 model at this same magnitude, our counts are consistent with either ~ 0.5 mag of luminosity 
evolution or a factor of ~ 1.5 enhancement in the galaxy density. While these models fail to reproduce 
the very faint Hubble Deep Field observations (see Ferguson & Babul (1998)), they do reasonably well at 



brighter magnitudes and thus provide a plausible comparison in the flux regime covered by our survey. 



4. The Two-point Correlation Function 



We use the |Landy k. Szalay (1993)| estimator, uj(9) = (DD—2DR+RR)/RR, to compute the two-point 



angular correlation function. Here DD is the observed number of galaxy pairs with separation between 9 
and 9 + 59, RR is the number of such pairs in a randomly distributed sample with identical boundaries 
and galaxy surface density, and DR is the number of such pairs in a cross-correlation between the observed 
catalog and a random realization. This is the minimum variance estimator and fully corrects for mask and 
edge effects up to scales approaching the survey size. Our results for 9 independent magnitude slices, each 
based on 100 random realizations, are shown in Figure || after the application of the integral constraint (IC) 
and stellar contamination corrections. The IC correction attempts to offset the artificial reduction in power 
on large scales that occurs when the mean galaxy density and uj{9) are determined from the same survey 



( Groth fc Peebles (1977) ). The small reduction in amplitude introduced by misclassified stars is corrected 
by multiplying the IC-corrected lu(9) values by the factor Nq^/ '(Nobj — Nstar) 2 where Nobj is the total 
number of objects in a given magnitude bin and N$tar is an estimate of the number of misclassified stars in 
the same bin. Since galaxies significantly outnumber stars when / > 21.5, stellar contamination effects are 
most significant for brighter magnitudes. However, that is precisely where the PSF-based classifier works 
extremely well. Hence we take the number of misclassified stars to be approximately equal to the square 
root of the number of stars in the bin. We have confirmed that the results from Landy-Szalay estimator are 
in excellent agreement with those from alternative estimators for uj(9) such as the counts-in-cells method 
or the ensemble estimator. In each magnitude slice, we only include data from those CCD images in the 
survey that are complete to at least 0.5 mag fainter than the upper mag limit of the slice. This assures that 
large-scale artificial power from frame-to-frame depth variations is eliminated. The depth of each frame is 
determined by identifying the location of the peak object counts (stars plus galaxies) and then subtracting 
0.25 mag. 

We find excellent agreement between our results for 17 < / < 18.25 and the APM lo{9) determination 
( Maddox et al. (1996)| ) derived from galaxies with 17 < bj < 20 (dashed line). Table 2 presents our cor- 



relation function data (multiplied by 100 for compactness) in each magnitude bin. The data in this table 
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have the IC and stellar contamination corrections applied. Table 3 presents the median / magnitudes, 
number of galaxies each subsample, best fit power law slopes, the corrected amplitudes at 0.5, 1, and 3 
arcminutes (interpolated from results in Table 2), and the values of the IC and stellar contamination correc- 
tions. The slope of the correlation function is determined from the best fit to the power law representation 
oj(&) = A W 9 S . The formal error in uj(9) is a u = + u;(#))/RR. This expression produces error values 
which are comparable to those estimated by the statistical bootstrap method. The errorbars in Figure || 
(and Table 2) include an often ignored additional error term, added in quadrature to the formal error, 
arising from uncertainty in the value of the IC. The IC uncertainty arises because one typically assumes a 
power-law form for ui(0) a priori and thus, on large scales, systematic errors can be introduced if this is not 
an accurate representation of the true correlation function. We estimate this uncertainty by measuring the 
change in the IC as the assumed power law slope is varied by ±0.1. Strictly, these two error components 
are not completely independent (and thus a quadrature addition is not completely proper). However, the 
quadrature addition provides a reasonable approximation and, for the present survey, the IC uncertainty 
is only important (i.e., becomes comparable to the formal error above) on scales 9 ^ 25 arcminutes. 



4.1. Dependence of the Amplitude on Magnitude 

Figure || shows the dependence of the amplitude of ui(9) at 1' scale on the median / magnitude of 



the subsample. Results from Efstathiou et al. (1991)|, |Campos et al. (1995) 


, Neuschaffer & Windhorst 


(1995), |Lidman & Peterson (1996) 


, Woods & Fahlman (1997) 


, and Brainerd & Smail (1997) 


are shown for 



comparison. We find a smooth decline in amplitude with median / magnitude that is consistent with the 
results from Campos et al. (1995) and Brainerd &: Smail (1997) . For I ,$ 21, our results are also consistent 



with those of [Lidman Sz Peterson (1996)| , although their results are systematically lower in amplitude by 



20% to 45%. At I > 22, we differ at more than the 3a level from the results of Lidman &: Peterson (1996) 



and |Woods fc Fahlman (1997) 



Our results are substantially more reliable than previous estimates over the range / < 24 because our 
large survey area overcomes the most significant problem faced by smaller surveys: the cosmic variance in 
galaxy surface density and clustering. For small surveys, these effects result in large fluctuations in the 
correlation amplitude and in a large value of the IC correction on arcminute scales. At large scales both of 
these effects depend on the average of the correlation function over the survey area. The ratio of variances 
expected between two different surveys is at least (9 1 /9 2 ) S/2 , where 6 X and 9 2 are the effective angular 
scales of the two surveys and S is the slope of the angular correlation function. While cosmic variance is 
boosted at small scales by the higher order correlations neglected in this simplified formula, the formula 
should be a reasonable accurate approximation for the ratio of finite volume errors. Thus, relative to many 
of the previous surveys, the cosmic variance on arcminute scales in our survey is suppressed by more than 
a factor of 3. Moreover, the continuity and simple geometry of our survey together with the large number 
of galaxies employed ensure that edge and discreteness effects are negligible over a large dynamic range 
of scales. The advantage of this large, contiguous area surveys is, therefore, even better than the above 
estimate. 

In addition to the above effects, which primarily depend on the variance only, the non-Gaussian error 
distribution for to(9) is likely to be skewed positively (Szapudi &: Colombi (1996)) as well. This can result 
in an apparent low bias in the correlation amplitude for smaller surveys since there are many areas of 
sky with values somewhat lower than the mean, which are not balanced by the relatively few areas with 
values much higher than the mean. We use our survey data to directly measure the true variance caused by 
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surface density and clustering properties in uj(1') in many independent but small contiguous areas, and find 
that variance is significant on 1 arcminute scales. We chose a cell size of 16' x 16' — corresponding to the 
field of view of a single KPNO 4m CCD exposure (which is still significantly larger than the independent 
fields of previous surveys). We analyze 250 such fields and computed uj(8) for 21 < I < 22, a magnitude 
range where we begin to see growing discrepancy between different surveys. The histogram of w(l') values 
is shown in Figure & We find that the scatter in u>(l') is comparable to its mean value (~ 0.045). Extreme 
values can reach 3x the mean. Thus, in order to measure w(l') with a S/N of 10 one requires at least 100 
such independent fields even at only moderately faint magnitudes. There is also a 20% offset between the 
full survey value and the mean of the 250 independent fields, but this difference is entirely consistent with 
uncertainties in the IC correction for the individual CCD frame results. 

The IC correction, like cosmic variance, is determined by the effective survey area, and to lesser 
extent, by the survey geometry. The largest contiguous area used in u>(9) measurements at I > 22 in the 
other works cited is 0.25 deg 2 (compared with 14.7 deg 2 here). The IC corrections at 1' in our survey 
are negligible ( ;$ O.O5u;(0)) whereas they can be a significant fraction (~ 50 — 100%) of the amplitude of 
oj(9) in some of the previous surveys (e.g., Efstathiou et al. (1991)). Hence, while the agreement between 



Prainerd fc Smail (1997) and our results is excellent, their two survey fields only cover about 0.02 deg 2 and, 



therefore, additional imaging over larger areas at fainter limits (/ > 24) is still desired. The |Campos et al 



(1995)Qesults are based on a survey of 64 11' x 11' non-contiguous fields (2.2 deg 2 total area) distributed 
over an 8° x 10° region. 



4.1.1. Models 



The dependence of the angular correlation function amplitude on magnitude can be modeled using 
the equation ( [Limber (1953] , Peebles (1980)1) 



u(0) 



7T- 



r[( 7 -l)/2] A(e) 



0T 



-1 o 



(1) 



r(7/2) 

that relates uj(8) to the spatial correlation function, £(r) = (r/r )~~ ( . A power-law spatial correlation 
function implies that it's two dimensional counterpart, oj(6), will also take on a power-law form. The 
function A(e) depends on the the galaxy redshift distribution, N(z), and on the evolution of £(r). If this 
evolution is parameterized as (e.g., Efstathiou et al. (1991) and [Woods fc Fahlman (1997) ) 



£ > (r,z) = (-)^(l + z) 
r n 



-(3+6) 



then A(e) is 



A(e) 



D(z)^g(z)-\l + Z )-^ 



dz 



fdN(z] 
Jo V dz 



dz 



-2 



where D(z), the angular diameter distance, and g(z) are 



c / q Q z + (q Q - l)( v / l + 2q Q z - 1) 



ql(l + zY 



g(z)- 1 = _£((! + z) 2 ^l + 2q z 



(2) 



(3) 



(4) 



(5) 



In this parameterization, two special cases are noteworthy. Clustering that remains constant in co-moving 
coordinates yields e = 7 — 3 ~ —1.2. Clustering that remains fixed with respect to physical coordinates 
yields e = 0. 
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We use several different redshift distributions in our fitting procedure (all shown in Figure ^). The 
first is based on fits of the empirical function dN(z)/dz oc z 2 exp(— (z/z ) 2 ) to the CFRS data and its 
extrapolation to fainter magnitudes (Lilly et al. (199"5)|; hereafter CFRS model). This parameterization 



provides a very good representation to the observed data ( [Brainerd fc Smail (1997)D and allows one to use 
a smoothly varying function in the Limber equation, thus removing shot noise effects present in the actual 
data. The median redshifts for the CFRS-based N(z) model are 0.095, 0.145, 0.260, 0.420, 0.540, 0.615, 
0.690, 0.755, 0.810, and 0.865 for I median mags of 16.5, 17.5, 18.5, 19.5, 20.5, 21.25, 21.75, 22.25, 22.75, 
and 23.25, respectively. The remaining N(z) models are generated by integrating an evolving Schechter 
luminosity function in successive redshift shells (hereafter EvLF models) and within magnitude limits that 
correspond to our chosen subsamples. The evolution is assumed to be a pure luminosity evolution as 
represented by the Schechter LF parameters a = —1.1, Mj(z) = —20.9 + 51og/i — (3z, and (3 = 1, 1.5,2. 
The (3 values were chosen to approximately bracket the CFRS model. All the models contain a significant 
fraction of high-z galaxies by I = 23.5: ~ 30 — 65% with z > 1 and ~ 3 — 30% with z > 1.5 for the EvLF 
models; 44% with z > 1 and 11% with z > 1.5 for the CFRS model. 

For each N(z) model, we find the best-fit to our u(l') vs I magnitude data by identifying the com- 
bination of r Q and e that minimize the x 2 statistic. Table 4 and Figures [/] and |8| summarize the re- 
sults of the fitting procedure. The best fit in the range 16 < I < 23 using the CFRS N(z) model gives 
r Q = 5.6/j- 1 Mpc and e = -0.20. The best fits for the EvLF N(z) models are (r , e) = (3.8k- 1 Mpc, -0.80), 
(3.8/1" 1 Mpc, -1.40), and (A.lh- 1 Mpc, -1.70) for (3 =1, 1.5, and 2, respectively. If we fix the correlation 



length to be 5.5 ± 1.5h 1 Mpc, a value that is typical of that found from local redshift surveys (e.g., de 



Lappar jmtTGeller, fc Huchra (19881 , [Jing, Mo, fc Boerner (1998)| , |Tucker et al. (1997)|) , the best fit e 



values fall in the range —0.4 < e < +1.3. Restricting the fits to data with I < 20, in general, also yields 
more positive e values and higher correlation lengths that are consistent with those from the above local 
redshift surveys (see columns 3 and 5 in Table 4). Fits to data in subsets between / = 18 and / = 23 are 
consistent with the results obtained using all data with / < 23. In other words, we see little change in 
the correlation length with magnitude for 18 < / < 23 for a given redshift distribution. All fits were done 
assuming a slope of 7 = 1.8 and q Q = 0.1 (for q Q = 0.5 the best fit r a values are ~ 7% lower). 

Figure || shows the reduced x 2 contours for the four N(z) models as functions of r Q and e. Two things 
are clear from the x 2 measurements. First, there is a significant covariance between r a and e - for a given 
redshift distribution, lower correlation lengths are coupled with lower e values. Second, the best-fit is 
located within a fairly shallow minimum at least in the direction of the r Q — e covariance. Thus, although 
the best fits prefer e < 0, positive values are included in the la contours, albeit with correspondingly 
higher correlation lengths. We provide an approximate description of the covariance between r Q and e 
by determining the best-fit line to the data within the la contour. The fit is determined by computing 
the mean between the fits r Q = /(e) and e = f(r ). The typical la uncertainty in r a at a given e is 
±0.15h- 1 Mpc and ±0.23/i _1 Mpc for the EvLF and CFRS models, respectively (for the full / < 23 
dataset). The typical la uncertainty in e at a given r Q is ±0.18 (for all models). These fit parameters are 
given in Table 4. 



4.2. Dependence of e on Angular Scale 



Our large survey allows us to perform model fits at scales other than 1'. We perform a similar analyses 
for 9 = 0.5', 3', 10', and 30'. The results are shown in Figure ^. Figure 1C shows the best fit e value as 
a function of angular scale obtained when the CFRS redshift model is used. The results for the EvLF 
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redshift models are similar shape although there is a zeropoint shift in the e values. For / > 22, there is 
a trend towards larger e values as the angular scale decreases, perhaps suggesting that the clustering on 
the smallest scales evolves more rapidly. The uncertainties in e on scales greater than 10 arcminutes are 
substantial, however, and the trend is not seen at I < 21. We show in the following section that there is a 
significant flattening of the slope of the correlation function when I > 22 as well. These two trends are not 
easily identified with any known observational or instrumental effects (see discussion below) but clearly a 
deeper, wide area survey is needed to confirm these results. 



4.3. Dependence of the Slope on Magnitude 

We find no significant dependence of the correlation function slope on magnitude for I < 22 as shown 



in Figure 11. Over this magnitude range and for 1' < < 20' the best fit slope is 5 = —0.80 ±0.02. For our 
two faintest bins (/ > 22), a flattening of the slope on angular scales 5' is seen with 5 = — 0.48±0.04 for 
22 < / < 23 (see Figure || and Table 3 as well). The trend is not explained by seeing or deblending effects as 
it is also seen in the subset of the data with the best resolution (FWHM < 1.25"). It could, in principle, be 
introduced by poorly calibrated frame to frame photometric zeropoint variations. However, our observing 
strategy enables us to measure and correct for these variations quite well (see §2) and any residual variations 
are not sufficient to introduce a flattening of this magnitude. Furthermore, the 22 < / < 22.5 bin is still 
1 magnitude brighter than the completeness limit. Experiments with simulated CCD data show that our 
object detection software is not dependent on the clustering properties of the objects. As an additional 
check, however, we compute lu(9) independently for the 4 quadrants of the survey, for the single run with 
the best atmospheric transparency, and for each individual CCD frame. We see the flattening in these 
subsamples as well. One would not expect an effect introduced by zeropoint miscalibration to survive 
all these experiments. This dependence has been seen in two other surveys at approximately the same 



magnitude and over the same angular scales as seen here. Neuschaffer fc Windhorst (1995) report slopes 
of 6 = —0.5 by g = 25 based on two independent fields, each one covering 0.25 contiguous deg 2 . Campog 



et al. (1995) find best fit slopes of 5 = —0.6 and 5 = —0.55 for their R < 22 and R < 23 samples, 
respectively. There is a physical model that can be invoked to explain the flattening; we will summarize it 
in the next section. While it is intriguing that our results are roughly consistent with those seen in these 
two independent surveys done at different wavelengths, we remain cautious for now since the effect is only 
seen in our faintest bins. Other large, contiguous surveys now underway (pey fc Jannuzi 1997, Falco 1997) 



should provide an important check on this result. We note that the constraints on r Q and e discussed above 
do not change if we exclude the two faintest points from the analysis. 



Discussion 



Our constraints on to(0) reveal a dependence of its amplitude on median I magnitude that appears 
to fit the model of the redshift evolution of £(r) parameterized as £(r, z) = (^) _7 (1 + z)~( 3+e \ In fact, 
this model provides an acceptable representation of the data over the range 16 < / < 25 if the results 



of Brainerd &: Smail (1997) are also included in the fits. The decline in the amplitude of the correlation 
function (over a wide range of scales) with / magnitude is about 3 times steeper over the range 16 ^, I ^ 19 
than over the range 20 ;$ / ;$ 23. The N(z) distributions that yield good fits to magnitude-u;(#) relation 



typically have about 40% of the galaxies at z > 1 by / = 23, consistent with the conclusions of |Brainerd 



fc Smail (1997) . A shallower decline in the amplitude of ui{9) with magnitude for / ^ 21 can be induced 
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by the reduction of the proper volume element at higher z. Brainerd fc Smail (1997) use this argument to 
support the assertion that a significant fraction of faint galaxies are at z £ 1- However, a high-z cutoff can 
also induce similar behavior since in this case the effective depth and redshift distribution are not changing 



dramatically as one goes to fainter limits. If such a cutoff were introduced, as noted by Lidman & Peterson 
(1996)[]by a selection bias against the detection of high-z galaxies a shallow decline in correlation length 
with magnitude would result. However, a strong redshift cutoff would also effect the number counts and the 



corresponding turnover in the counts is not seen (see also Woods &: Fahlman (1997)1 ), Given the excellent 



agreement we find with the HDF number counts, we conclude that a selection bias of this nature is not a 
significant problem for our sample. 

Taken at face value, our best fits to the above model suggest the evolution of the spatial correlation 
function is relatively mild over the depth of the survey. When we allow r to be a free parameter in the 
fit, we obtain exponents in the range —1.7 ^ e ^ — 0.2 with 3.8 ^ r Q ^ 5.6/i _1 Mpc depending on the 
assumed redshift distribution. Growth in the clustering of galaxies (in proper coordinates) requires e > 0. 
The 95% confidence limits are substantial, nevertheless, given the strong covariance between r Q and e. As 
Figure || demonstrates, our data are also consistent with e > providing 4.5 iS r a ;$ 7/i~ 1 Mpc. If e ^ 1 
then values of r a 4/i -1 Mpc are strongly rejected. As already noted, redshift surveys of local galaxies 
find r a ~ 5.5 ± 1.5/i -1 Mpc. 

The interpretation of the goodness of the fit to the above model, however, must take into account a 
number of additional selection effects that may mask more substantial evolution of the correlation function. 
It is known that low-z galaxy clustering depends on both the morphology and luminosity of the objects 
being studied (e.g., |Davis fc Geller (1976)| , Moore et al. (1994)| , |Loveday et al. (1995)| , |Guzzo et al. (1997)1 , 



Valotto fc Lambas (19971 ). The sense of these trends is that elliptical galaxies and more luminous galaxies 



tend to have larger correlation lengths than spiral galaxies and less luminous galaxies, respectively. We also 
know that the mean absolute luminosity of the objects in a flux-limited sample will increase with redshift. 
The relative insensitivity of our derived r Q values to I magnitude for galaxies with / > 17 could therefore 
be, in part, the result of two competing effects. As one goes fainter the mean redshift of the galaxies in 
the survey tends to increase. If galaxies at 0.5 £ z 1 exhibit a similar luminosity dependent clustering 
as local galaxies and/or /—band selection enhances the elliptical galaxy fraction at these redshifts then 
it is possible that a decrease in the spatial correlation length with magnitude could be partially offset by 
the tendency for more luminous and/or early type galaxies to be more strongly clustered. This could also 
explain the observed preference for negative e values (which suggests only modest clustering evolution). 
The precise evolution of e and r Q with redshift is intimately connected to the density pertubation power 
spectrum and the galaxy merger rate ( Moscardini et al. (1997)] ). The complexity of the dependence makes 



a unique theoretical interpretation the observational constraints problematic. None the less, the constraints 
are now solidly determined on degree-scales and less down to I = 23 from the present survey. 

More stringent constraints on clustering evolution will require the addition of multiple passbands or 
redshift data. Efstathiou et al. (1991) report a marginal increase in the clustering amplitude as one selects 



redder passbands and posit that normal galaxies dominate the composition of fainter /—band selected 
galaxy samples as opposed to the weakly clustered faint galaxies, that appear to dominate U and B 
selected surveys. Lidman Sz Peterson (1996)] also find a marginally significant color dependence, showing 



that red (V — I > 1.5) galaxies have a larger clustering amplitude than blue galaxies. |Le Fevre et al 
(1996)Qind no color dependence to the clustering properties of galaxies at z > 0.5 from an analysis of the 
CFRS data. They have also measured the spatial correlation function directly from their redshift data. 
They find a correlation length of r a = 1.57 ± 0.09/i -1 Mpc at z ~ 0.5 which, if they are sampling a mix 
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of galaxies similar to those in local redshift surveys, implies < e < 2. Our results are consistent with 
this e range. From equation (2), we know r (z) = r o (0) (1 + z) ( - (3+e)/7) . Using our best-fit results for the 
CFRS redshift distribution we, thus, find that r(z = 0.5) « 3.0/1" 1 Mpc. Results for the EvLF redshift 
distribution models yield r(z = 0.5) = 2.7 ± 0.4/i _1 Mpc. The CFRS galaxies were also /—band selected 
so any difference in results between the CFRS and our angular survey are due either to projection effects 
(which tend to wash out clustering) or the volume sampled as opposed to intrinsic differences between the 
galaxy population. 

I I 

Connolly, Szalay, & Brunner (1998) determined the redshift dependence of the amplitude of oj{0) 

using photometric redshifts derived from the HDF survey. For z > 0.4, they find r Q = 2.37/i _1 Mpc 
and e = — 0.41q'q5 However, in order to be consistent with low-z surveys, they would require e ~ 2. 
They suggest, therefore, that the expression for the evolution of the spatial correlation function, £(r, z) = 
(^-) -7 (1+ z)~^ +€ \ is not particularly good. Others have made similar claims based on poor fits to the oj{6) 
vs. magnitude relation. This is in contrast with the excellent fit we obtain extending over 7 magnitudes 
(Figure ^) . We speculate that the discrepancy is caused by cosmic variance due to the small area used in 
previous measurements of the correlation function. Connolly, Szalay, & Brunner (1998) note that the HDF 



survey subtends only 800/i _1 kpc at z = 1 (q D = 0.1), a fraction of the galaxy correlation length. The CFRS 
is larger (five 10' x 10' fields) but still covers a relatively small volume (10' spans 3.1/j^ 1 Mpc at z = 1). 
Hence, these surveys are also subject to the same sorts of problems associated with cosmic scatter and 
undersampling of LSS as small area angular surveys. Indeed, Be Lapparent, Geller, &: Huchra (1988)| make 
a clear demonstration of this point at low-z: even in their redshift survey of ~ 1810 galaxies (z < 0.05) 
filling a volume of 5 x 10 5 /i -3 Mpc 3 , the density fluctuations caused by LSS prevent the determination of r a 
by better than a factor of 2. Given that the deep, redshift probes sample roughly comparable volumes (but 
with the added complication that they stretch over a significant fraction of cosmic time) , it seems optimistic 
to expect reliable constraints to be determined from the existing surveys. In fact, it is quite remarkable 
that the values obtained from them are in agreement with the local correlation length to within a factor 
of two. Not until much larger volumes are surveyed will the constraints on r Q at high redshift be robustly 
measured. If the results from direct measurements of the spatial correlation function at higher redshifts 
are ultimately accurate, however, and assuming their sample composition is not significantly different from 
the low redshift surveys, then the basic model may be in need of refinement. 



One possible refinement, discussed in some detail already by Neuschaffer fc Windhorst (1995) , is to 
include a redshift-dependent spatial correlation function slope. This would be a natural consequence of 
models in which there is a scale dependence to the growth of structure. For example, if galactic scale 
structures grow significantly faster than structures on scales ^ 20h~ 1 Mpc the correlation function slope 
would steepen with time. This would also cause the slope of to(8) to flatten at fainter magnitudes on 
small scales and, depending on the details of the density perturbation power spectrum and the merger 
rate, could also result in an increase in e with redshift on small scales. Adopting the parameterization 



j(z) = 1.8(1 + z me d) c proposed by Neuschaffer fc Windhorst (1995)| , implies C = 0.35 ± 0.10 which is 



consistent with their limit. As |Neuschaffer fe Windhorst (1995)| demonstrate, values in the range C ^ 0.4 



are consistent with structure formation models which have significant power on co-moving scales of ^ 10ft, -1 
Mpc but disfavor models with little small scale power such as HDM. However, since the flattening is only 
seen at / ^ 22, the above parameterization of the slope redshift dependence is not particularly good - there 
appears to be no real redshift dependence to the slope until z mei i ^ 0.6. 
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6. Summary 

We have produced an /—band selected catalog of ~ 710, 000 galaxies for studying the evolution of 
structure out to z ~ 1. Our survey puts new and stringent limits on the magnitude dependence of the 
two-point angular correlation function. The survey is the first of a new generation of deep, wide-area 
imaging surveys made feasible by large-format CCDs and mosaic cameras. The key results presented in 
this paper are 



1. Galaxy number counts are consistent with modest evolution in luminosity (0.5 — 1 mag, depending 
on Q Q ) and/or density (factor of 1.5 — 2). The counts are consistent with the F814W HDF counts in 
the range 21 < / < 24 suggesting that the HST-based catalog does not contain a substantial number 
of sub-galactic components in this flux range. 

2. The two-point correlation function agrees remarkably well with that of the APM in the magnitude 
range common to both surveys. At the faint end, our measurements of oj{9) are consistent with those 



of Brainerd &: Smail (1997) , who have measured u)(9) down to / = 25. The large contiguous size of 
our survey minimizes the effects of cosmic scatter and provides a firm determination of ui(9) over 7 
magnitudes. 

3. The amplitude of the two-point angular correlation function decreases monotonically with increasing 
magnitude over the range 16 < / < 23. The decline is well fit by a redshift dependence of the spatial 
correlation function parameterized as (1 + z) _ ( 3+e ). Our best fit correlation lengths for galaxies with 
I < 20 are consistent with those from low-z redshift surveys. At fainter magnitudes we find correlation 
lengths in the range 3.8 < r a < 5.5/i _1 Mpc depending on the assumed redshift distribution. If e ^ 1 
then we strongly reject r Q 4/i _1 Mpc for our / > 20 sample. However, more negative e values are 
preferred at I > 20, suggesting the effective clustering is not evolving as rapidly as linear perturbation 
theory predicts. The possible preference for inclusion of more luminous and/or early type galaxies 
at fainter magnitudes in the survey could explain this trend. 

4. Our derived correlation length at z ~ 0.5 is about a factor of 2 larger than that derived from the 
CFRS. While a complete physical interpretation of our results requires additional color and redshift 
data, it is also clear that the volumes sampled by existing deep, redshift probes are inadequate to 
properly constrain the two-point spatial correlation function to better than the above difference. If 
the galaxy population sampled at I ^ 20 is similar to that in the CFRS, then our deprojection of 
£(r) may indeed be more accurate by virtue our reduced sensitivity to cosmic scatter. 

5. For / < 22, the mean value of power-law slope of uj(0) is — 0.80±0.02 and is independent of magnitude. 
We find a slope of 5 ~ —0.5 is a better fit to the I > 22 data, however. This flatter slope is consistent 



with the results of Campos et al. (1995) and Neuschaffer &; Windhorst (1995) . We also detect an 



increase in e on scales less than 1 arcminute for / > 22. We can identify no instrumental or software- 
related cause for these effects but remain cautious in interpreting the full signals as real phenomena. 
Such trends are, however, consistent with structure formation models in which small scale power 
increases more rapidly than power on scales larger than ~ 10/i -1 Mpc. 
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Fig. 1. — The differential galaxy surface density as a function of isophotal / magnitude. The /—band 
counts from the Hubble Deep Field (open circles; Williams et al. (1996) ) are shown as well. The models 
(dashed and solid lines) are from [Ferguson fc Babul (1998) . 
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Fig. 2. — The differential galaxy surface density as a function of isophotal I magnitude for various surveys. 
Counts have been normalized by dex(0.6(7 — 18)) to enhance small variations and departures from the 
expectations of Euclidean space. 
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Fig. 3. — The two-point angular correlation functions, normalized by #~ ' 75 , for 17 < I < 23. Results for 
each magnitude range are vertically offset by 0.5 dex, as well, for clarity (the value of each offset is shown 
explicitly on plot). The best-fit power law for the APM 17 < bj < 20 galaxy sample is shown as a dashed 
line (its vertical offset is the same as that for the 17 < I < 18 sample). An arrow denotes the CCD field 
of view. 
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Fig. 4. — Amplitude of uj(6) at 1' vs. median / magnitude. Our results are the filled circles. Other recent 
measurements are also shown. 
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Fig. 5. — The histogram of ui(V) in the range 21 < / < 22 for 250 independent 16' x 16' CCD fields. The 
mean value is 0.045 but the variance is substantial owing to significant cosmic scatter on this scale at these 
magnitudes. The ~ 20% offset between the mean u>(l') for the independent CCD fields and that for full 
survey is entirely consistent with the uncertainties in the small field IC correction. 
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Fig. 6. — The differential redshift distributions for the CFRS model and the evolving luminosity function 
(EvLF) models for the magnitude intervals 16 < / < 17, 17 < I < 18, 18 < / < 19, 19 < / < 20, 
20 < / < 21, 21 < I < 22, and 22.5 < I < 23.5. The EvLF models assume a Schechter luminosity function 
with a = -1.1, M*j{z) = -20.9 + 51og/i - 0z, and (3 = 1, 1.5, 2. 
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Fig. 7. — Best fits to our u>(l') vs / mag data using a Limber deprojection and various plausible redshift 
distributions. In the case of the EvLF (3 = 1.5 N(z) model, we show fits to all the data and to just data 
with / < 20. Fits to data with / > 20 give similar results as those done using all data. 




Fig. 8. — The reduced x 2 contours for the fits to equation 1 using the EvLF (f3 = 1,1.5,2) and CFRS 
redshift distributions. Results shown here are for correlation amplitudes at 1 arcminute and I me dian < 23. 
The la (heavy line), 2a, and 3a contours are shown. 




Fig. 9. — The amplitude of lj(0r) as a function of median / magnitude for Or = 0.5', 1', 3', 10', and 30'. 
Data have been corrected for integral constraint and stellar contamination. Best-fit models using the EvLF 
P = 1.5 and CFRS redshift distributions are shown as solid and dashed lines, respectively. 
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Fig. 10. — The best-fit e value as a function of angular scale for data with I < 21, I < 22, and I < 23. 
The results shown are for the CFRS redshift distribution model. Results for the EvLF models are similar 
in shape but are offset vertically relative to the CFRS fits. 
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Fig. 11. — The best fit power law slope of the two-point correlation function as a function of / magnitude. 
The results for two different angular fitting limits are shown. The dashed line shows the mean slope derived 
for I < 22. 



Table 1. Differential I-band Galaxy Counts 



I mag log(N) a log(N Co r) a Vhigh viow h I mag log(N) a log(N Cor ) a @high &low 



13.25 


-0 


5671 


-0 


5671 





3010 





5671 


18 


75 


2 


9433 


2 


9433 





0076 





0077 


13.50 


-0 


2661 


-0 


2661 





2323 





5333 


19 


00 


3 


0641 


3 


0641 





0066 





0067 


13.75 





3360 





3360 





1315 





1895 


19 


25 


3 


1613 


3 


1613 





0059 





0060 


14.00 





3872 





3872 





1249 





1761 


19 


50 


3 


2696 


3 


2696 





0052 





0053 


14.25 





4329 





4329 





1193 





1651 


19 


75 


3 


3633 


3 


3633 





0047 





0047 


14.50 





7551 





7551 





0857 





1069 


20 


00 


3 


4547 


3 


4547 





0042 





0043 


14.75 





8801 





8801 





0752 





0910 


20 


25 


3 


5482 


3 


5482 





0038 





0038 


15.00 


1 


0240 


1 


0240 





0645 





0758 


20 


50 


3 


6360 


3 


6360 





0034 





0035 


15.25 


1 


1142 


1 


1142 





0586 





0677 


20 


75 


3 


7218 


3 


7218 





0031 





0031 


15.50 


1 


2590 


1 


2590 





0501 





0566 


21 


00 


3 


8007 


3 


8007 





0028 





0029 


15.75 


1 


3967 


1 


3967 





0431 





0478 


21 


25 


3 


8904 


3 


8904 





0026 





0026 


16.00 


1 


5665 


1 


5665 





0357 





0389 


21 


50 


3 


9690 


3 


9690 





0023 





0023 


16.25 


1 


6582 


1 


6582 





0323 





0349 


21 


75 


4 


0486 


4 


0486 





0021 





0021 


16.50 


1 


7870 


1 


7870 





0280 





0299 


22 


00 


4 


1196 


4 


1437 





0019 





0019 


16.75 


2 


0046 


2 


0046 





0219 





0231 


22 


25 


4 


1825 


4 


2000 





0027 





0027 


17.00 


2 


0985 


2 


0985 





0197 





0207 


22 


50 


4 


2579 


4 


2820 





0024 





0024 


17.25 


2 


2175 


2 


2175 





0173 





0180 


22 


75 


4 


3222 


4 


3554 





0022 





0022 


17.50 


2 


3327 


2 


3327 





0151 





0157 


23 


00 


4 


3873 


4 


4385 





0020 





0020 


17.75 


2 


4875 


2 


4875 





0127 





0131 


23 


25 


4 


4549 


4 


5438 





0018 





0018 


18.00 


2 


5600 


2 


5600 





0117 





0120 


23 


50 


4 


5166 


4 


6383 





0016 





0016 


18.25 


2 


7164 


2 


7164 





0098 





0100 


23 


75 


4 


5768 


4 


7282 





0015 





0015 


18.50 


2 


8314 


2 


8314 





0086 





0088 























a N and Ncor are the differential galaxy counts before and after correction for survey incomplete- 
ness, respectively. The units are galaxies deg -2 mag -1 . 

b The 1-sigma Poisson errors in log(Ncjor)- 



Table 2. Angular Two- Point Correlation Function 



e 


I < 17 


17 < I < 18 


18 < J < 19 


19 < / < 20 


20 < / < 21 


21 < / < 21.5 


21.5 < / < 22 


22 < / < 22.5 


22.5 < / < 23 


0.37' 


192 ±84.5 


136 ±33.8 


57.4 ±8.96 


29.4 ±3.02 


13.1 ± 1.23 


8.56 ± 1.36 


6.46 ±0.94 


6.87 ±0.98 


7.28 ±0.72 


0.47 


79.2 ±53.9 


89.2 ±23.4 


48.3 ±6.94 


21.7 ±2.32 


12.2 ± 0.98 


10.3 ± 1.10 


9.44 ±0.77 


4.08 ±0.77 


3.96 ±0.57 


0.59 


131 ± 46.8 


74.0 ± 17.8 


23.7 ±5.00 


19.2 ± 1.85 


10.7 ±0.78 


7.86 ± 0.86 


7.16 ±0.61 


4.31 ±0.62 


4.14 ±0.46 


0.74 


167 ±40.2 


59.2 ± 14.1 


27.2 ±4.08 


15.1 ± 1.44 


10.6 ± 0.62 


5.96 ±0.68 


5.90 ±0.48 


4.60 ±0.50 


4.39 ±0.37 


0.94 


146 ±29.7 


63.5 ± 11.3 


21.2 ±3.19 


16.2 ± 1.17 


8.96 ± 0.50 


6.16 ±0.55 


5.64 ±0.39 


4.90 ±0.41 


4.34 ±0.30 


1.18 


115 ±22.4 


34.3 ±8.20 


16.5 ±2.51 


10.7 ±0.92 


7.18 ± 0.40 


4.63 ± 0.44 


4.61 ±0.31 


3.43 ±0.33 


4.04 ±0.25 


1.48 


95.7 ± 16.9 


44.7 ±6.76 


16.1 ±2.00 


10.4 ±0.74 


5.85 ± 0.33 


4.97 ±0.35 


3.92 ±0.25 


3.99 ±0.27 


3.23 ±0.21 


1.87 


84.1 ±13.3 


30.6 ± 5.16 


13.4 ± 1.59 


8.74 ±0.60 


5.19 ± 0.28 


4.42 ±0.29 


3.52 ±0.21 


3.07 ± 0.23 


2.99 ±0.17 


2.35 


67.1 ± 10.2 


29.0 ±4.12 


11.2 ± 1.27 


6.70 ±0.50 


4.17 ±0.23 


2.76 ±0.23 


2.96 ±0.17 


2.83 ±0.19 


2.70 ±0.15 


2.96 


47.7 ± 7.82 


21.6 ±3.23 


11.2 ± 1.03 


5.59 ± 0.42 


3.91 ± 0.20 


2.75 ±0.19 


2.51 ±0.15 


2.41 ±0.17 


2.48 ±0.13 


3.72 


54.8 ±6.50 


26.1 ±2.66 


9.51 ±0.84 


4.59 ± 0.36 


3.24 ±0.18 


2.47 ±0.16 


2.03 ±0.13 


2.39 ±0.15 


2.09 ±0.12 


4.69 


44.1 ± 5.17 


17.9 ±2.11 


7.72 ±0.69 


3.71 ± 0.31 


2.20 ±0.17 


2.12 ±0.14 


1.80 ±0.11 


1.92 ±0.14 


1.76 ±0.11 


5.90 


27.1 ±4.06 


17.6 ± 1.75 


5.17 ±0.58 


3.25 ± 0.28 


2.09 ± 0.16 


1.53 ±0.13 


1.77 ±0.10 


1.53 ±0.13 


1.31 ±0.11 


7.43 


30.2 ± 3.49 


12.7 ± 1.45 


4.89 ±0.50 


2.64 ± 0.26 


1.66 ± 0.15 


1.47 ±0.11 


1.28 ±0.10 


1.44 ±0.12 


1.21 ±0.10 


9.35 


17.7 ±2.93 


8.81 ± 1.23 


3.64 ± 0.43 


2.22 ± 0.25 


1.43 ± 0.14 


1.20 ±0.11 


0.98 ±0.09 


1.23 ±0.12 


0.94 ±0.10 


11.77 


16.9 ±2.61 


6.92 ± 1.08 


2.94 ±0.39 


1.95 ± 0.24 


1.28 ± 0.14 


0.84 ±0.10 


0.75 ±0.09 


0.92 ±0.12 


0.71 ±0.10 


14.82 


10.3 ±2.36 


5.22 ±0.97 


2.92 ±0.36 


1.53 ± 0.23 


0.94 ±0.14 


0.60 ±0.10 


0.59 ±0.09 


0.69 ±0.11 


0.52 ±0.10 


18.66 


10.0 ±2.20 


5.54 ±0.90 


1.81 ±0.34 


1.27 ±0.23 


0.75 ± 0.14 


0.44 ± 0.09 


0.40 ±0.09 


0.52 ±0.11 


0.36 ±0.10 


23.49 


11.1 ±2.10 


5.72 ±0.86 


1.61 ±0.33 


0.95 ± 0.22 


0.64 ±0.14 


0.45 ± 0.09 


0.37 ±0.09 


0.53 ±0.11 


0.27 ±0.10 


29.57 


7.79 ± 2.03 


3.87 ±0.82 


0.91 ±0.32 


0.63 ± 0.22 


0.38 ± 0.14 


0.27 ±0.09 


0.26 ±0.09 


0.51 ±0.11 


0.22 ±0.10 


37.22 


7.51 ± 1.99 


2.60 ±0.80 


1.01 ±0.32 


0.57 ±0.22 


0.41 ± 0.14 


0.23 ±0.09 


0.14 ±0.09 


0.35 ±0.11 


0.21 ±0.10 


46.86 


5.49 ± 1.96 


2.88 ±0.79 


0.91 ±0.31 


0.69 ± 0.22 


0.35 ± 0.14 


0.12 ±0.09 


0.05 ±0.09 


0.13 ±0.11 


0.25 ±0.10 


58.99 


2.21 ± 1.94 


2.37 ±0.78 


0.91 ±0.31 


0.67 ±0.22 


0.33 ±0.14 


0.12 ±0.09 


-0.07 ±0.09 


0.23 ±0.11 


0.05 ±0.10 


74.27 


4.96 ± 1.92 


1.75 ±0.77 


0.67 ±0.31 


0.60 ±0.22 


0.23 ±0.14 


0.11 ±0.09 


-0.06 ±0.09 


0.30 ±0.11 


0.06 ±0.10 


93.50 


5.94 ± 1.92 


0.79 ±0.77 


0.44 ±0.31 


0.53 ± 0.22 


0.24 ±0.14 


0.23 ± 0.09 


0.04 ±0.09 


0.40 ±0.11 


0.12 ±0.10 



a Divide u>(0) results by 100 to get actual values 



Table 3. Measured u(9) Parameters 



Imedian 


N gal 


<5 a 


5 h 


uj(0.5')cov C 


u{l')cor C 


Uj{3') C or C 


IC Corr. 


C s Corr. d 


16.46 


1443 


-0.879 ± 0.058 


-0.688 ± 0.087 


202 ± 101 


137 ±52.0 


56.4 ± 15.6 


0.0379 


1.167 


17.64 


3245 


-0.807 ±0.054 


-0.739 ± 0.079 


80.4 ±41.2 


53.9 ±19.5 


25.0 ± 6.46 


0.0151 


1.108 


18.62 


9987 


-0.796 ± 0.043 


-0.690 ±0.061 


30.1 ± 11.9 


22.5 ± 5.70 


9.86 ±2.06 


0.0061 


1.069 


19.61 


27246 


-0.809 ± 0.034 


-0.808 ± 0.041 


23.7 ±4.17 


13.6 ±2.09 


5.56 ±0.83 


0.0044 


1.046 


20.59 


62913 


-0.783 ±0.028 


-0.739 ± 0.030 


13.6 ± 1.75 


8.30 ±0.91 


3.55 ±0.41 


0.0027 


1.033 


21.27 


55245 


-0.788 ±0.033 


-0.639 ± 0.041 


8.91 ± 1.96 


6.19 ±0.99 


2.78 ± 0.39 


0.0018 


1.025 


21.77 


78082 


-0.742 ± 0.029 


-0.674 ± 0.034 


8.22 ± 1.37 


5.26 ± 0.70 


2.48 ± 0.39 


0.0017 


1.022 


22.27 


106975 


-0.627 ±0.034 


-0.475 ± 0.042 


5.89 ± 1.39 


4.55 ± 0.74 


2.44 ± 0.34 


0.0014 


1.020 


22.77 


137538 


-0.697 ±0.034 


-0.483 ± 0.035 


4.05 ± 1.03 


4.19 ±0.55 


2.25 ±0.27 


0.0016 


1.019 



a Slope for the best fit power law in the range 1' < 9 < 20'. 
b Slope for the best fit power law in the range 0.5' < 9 < 5'. 

c Divide by 100 to get actual amplitude. Based on a local interpolation to provide a more accurate estimate of the 
amplitude. The IC and stellar contamination (C s ) corrections have been applied. Corrected u>(9) = C s (oo(9) b s + IC) 

d Stellar contamination correction factor, C s = NQ b -/{Nobj — Nstar) 2 - 



Table 4. Best-fit r Q and e values (1 arcminute scale) 







r 


a To at 




a e a dX 




Covar. 


Covar. 






Sample 


N(z) Model 


(h- 1 Mpc) 


fixed e 


e 


fixed r 


2 b 


slope 


Intercept 


r Q Range d 


e Range d 


I < 20 


CFRS 


4.70 


0.32 


-1.50 


0.42 


0.11 


0.768 


5.822 


3 to 7 


-3 to +1 


/ < 20 


EvLF = 1 


5.30 


0.35 


+2.30 


0.64 


0.29 


0.555 


4.577 


3 to 7 


-3 to +4 


/ < 20 


EvLF =1.5 


5.20 


0.33 


+ 1.20 


0.66 


0.30 


0.495 


4.209 


3 to 7 


-3 to +4 


/ < 20 


EvLF = 2 


5.30 


0.36 


+0.40 


0.57 


0.30 


0.624 


5.060 


3 to 7 


-3 to +4 


/ < 23 


CFRS 


5.60 


0.23 


-0.20 


0.17 


0.24 


f.356 


5.874 


4 to 7 


-1.5 to +1 


/ < 23 


EvLF = 1 


3.80 


0.15 


-0.80 


0.2f 


0.39 


0.713 


4.389 


2.7 to 5 


-2.5 to +1 


/ < 23 


EvLF 0=1.5 


3.80 


0.15 


-1.40 


0.f8 


0.35 


0.830 


4.954 


2.8 to 5.1 


-2.7 to 


/ < 23 


EvLF = 2 


4.10 


0.16 


-1.70 


0+6 


0.26 


1.025 


5.798 


3 to 5.2 


-2.8 to -0.5 


/ < 23 


CFRS 


5.50 Fixed 




-0.30 


0.75 


0.25 










/ < 23 


EvLF = 1 


5.50 Fixed 




+ 1.30 


0.90 


2.09 










/ < 23 


EvLF 0=1.5 


5.50 Fixed 




+0.50 


0.80 


2.06 










/ < 23 


EvLF = 2 


5.50 Fixed 




-0.40 


0.70 


1.83 











Uncertainties in e at fixed r are the la values except when r is kept fixed. In that case, the errors shown are the 2a values. 
b The reduced x 2 value. 

c The slope and intercept for the best-fit line to the la x 2 contour. The line is r a — (slope)e + intercept. Values are the average 
parameters derived from the fits r a = /(e) and e = f(r a ). 

d The r a and e range over which the fit to the la \ 2 contour is valid. 



